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Abstract. We present a HI synthesis imaging study of NGC 864 (CIG 96), a spiral galaxy well isolated from similarly sized 
companions, yet presenting an intriguing asymmetry in its integral HI spectrum. The asymmetry in the HI profile is associated 
with a strong kinematical perturbation in the gaseous envelope of the galaxy, where at one side the decay of the rotation curve 
is faster than keplerian. We detect a small (M(HI) = 5 x 10 6 M G ) galaxy with a faint optical counterpart at ~ 80 kpc projected 
distance from NGC 864. This galaxy is probably not massive enough to have caused the perturbations in NGC 864. We discuss 
alternatives, such as the accretion of a gaseous companion at a radial velocity lower than the maximum. 
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1. Introduction 

The origin of asymmetries in isolated galaxies is not well un- 
derstood. Are they always induced by the presence of small 
companions? What is the influence of recent captures, and 
how long does it take for a parent galaxy to become again 
axisymmettic afterwards? Could some originate from internal, 
i as yet not well studied, long-lived dynamical instabilities (e.g. 
Baldwin et al. 1980)? Integral HI spectra provide a powerful 
statistical approach to this issue since they contain information 
about both the HI density distribution and the velocity field. 
We have observed, or compiled from the literature, single-dish 
21-cm line spectra for almost 800 galaxies (Espada et al. in 
prep.) as part of the AMIGA project (Analysis of the interstellar 
Medium of Isolated GAlaxies; Verdes-Montenegro et al. 2005), 
which aims to provide a baseline against which environmental 
effects can be evaluated. We have constructed a complete sam- 
ple of the most isolated galaxies, starting from the Catalogue of 
Isolated Galaxies (CIG, Karachentseva 1973; 1050 galaxies), 
selected on the basis of the distance to the nearest similarly 
sized galaxies. About 50% of these galaxies show asymmetric 
HI profiles, consistent with the results from previous work on 
smaller samples of isolated galaxies (104 galaxies in Haynes 
et al. 1998; 30 galaxies in Matthews et al. 1998). Surprisingly, 
other samples of galaxies in denser environments show only 
slightly larger values of this rate (50 - 80%, e.g. Swaters et al. 
2002, Richter & Sancisi 1994, Sulentic & Arp 1983). 
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HI synthesis imaging of very isolated galaxies with a highly 
lopsided HI profile may clarify the origin of these asymme- 
tries. Not only Hi-rich companions can be identified, but the 
presence of tidal features can be revealed, which would trace a 
past interaction. We have selected from our database a well de- 
fined sample of the most isolated galaxies showing significant 
asymmetties in their HI profiles (A„ >1.1, see Haynes et al. 
1998 for the definition), and mapped them with the VLA 1 in the 
21-cm HI line in order to look for signs of external interaction 
by analyzing in detail the HI distribution and its kinematics. 

One of the most interesting galaxies in our VLA sample 
is NGC 864, number 96 in the CIG catalog, which, as we 
shall see, shows a peculiar HI envelope in our 21-cm maps. 
It is classified as SAB(rs)c by de Vaucouleurs et al. (1991). It 
has an optical size of 4'. 7 x 3'5 and an apparent magnitude 
trig = 11.6. It has a very asymmetric HI profile (cf. Fig.Q. 
The systemic heliocentric velocity that we have measured for 
this galaxy (see Sec.|3j is 1561.6 km s _1 , corresponding to a 
distance of 17.3 Mpc (Ho = 75 km s _1 Mpc -1 ) after apply- 
ing the correction w.r.t. the CMB reference frame. This gives 
an optical luminosity of 1.1 x 10 10 L . The closest galaxy to 
NGC 864 at the same redshift is a small companion that we 
detected in our HI observations (see Sec.|2j with v/, e / = 1605 
km s _1 , located at a distance d ~ 80 kpc (153) of NGC 864, 
visible in the POSS2 as a faint object (identified in NED as 
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2MASX J02162657+0556038 with a reported magnitude of 
m B - 16.38) with an approximate size of 0'8 x 0'6 (physical 
diameter ~ 4 kpc). The other galaxies within a projected radius 
of d = 0.5 Mpc from NGC 864 are fainter by 3 to 7 magnitudes. 

2. HI synthesis observations and results 

Observations of the 21 cm line for NGC 864 were made with 
the VLA in its D configuration in July 2004 with 26 antennas. A 
bandwidth per IF of 3.125 MHz (from 1249.5 to 1895.2 km s" 1 ) 
was used, in 2IF correlator mode, giving a velocity resolution 
of 48.8 kHz (10.4 km s" 1 ) for the 64 individual channels after 
Hanning smoothing. The data were calibrated using the stan- 
dard VLA procedure in AIPS and were imaged using IMAGR. 
The average of the line-free channels has been subtracted from 
all the individual channels. The galaxy was detected between 
1436.7 km s _1 and 1686.6 km s _1 . The synthesized beam was 
49"8 x 46"2 (ax5). The rms noise level achieved after 4 hours 
is 0.66 mJy/beam. Primary beam correction has been applied 
to our maps. In Fig.@].a) and Fig.@].b) we display the channel 
maps smoothed to a beam size of 70"4 x 65"3, containing the 
HI emission at the observed radial velocities, from 1426.3 to 
1707.4 km s . We detect the small companion in the HI map 
from 1572.0 to 1655.3 km s" 1 at ff(2000.0) = 02 /, 16"'26: , 9 and 
5(2000.0) = 05°56'24:'0 . 

The total spectrum has been obtained by integrating the 
emission in the individual channel maps (solid line in Fig. ^ 
top). We have measured a HI content of M H i = 7.53 x 10 9 M G , 
in very good agreement with the single dish profile obtained 
by Haynes et al. (1998, dotted line in Fig.^ top), showing that 
there is no loss of flux in the synthesis imaging. The profile of 
the small companion is shown at the bottom of the same fig- 
ure. Subtraction of the flux of this dwarf galaxy from the total 
spectrum will not change its (asymmetric) shape. 

In Fig. 13 we show the integrated HI column density dis- 
tribution of NGC 864 (left) and the HI velocity field (right), 
both overlaid on a POSS2 red band image, the deepest one 
we found. These maps have been calculated as follows. The 
channel maps were smoothed with a gaussian tapering func- 
tion, and then non-signal pixels of the maps were blanked out. 
The channel maps where emission was detected were added up 
to produce the integrated emission map. The intensity weighted 
mean radial velocity field has been obtained in a similar way. 
The small companion mentioned above is clearly seen, together 
with its faint optical counterpart. 

The integrated emission map shows an unresolved HI de- 
pression in the center of the galaxy, while the highest column 
densities are located in a pseudoring structure with a size of ~ 
90" x 55", better seen in the grey scale map shown in Fig. [3] 
This seems to trace the outer spiral structure of NGC 864, al- 
though the spatial resolution of the HI data is too small for 
a detailed comparison. An arm-like feature seems to join this 
structure at a second enhancement in the column density dis- 
tribution. This enhancement has a ring-like shape with an ap- 
proximate size of 8'2 x 4'6 (Fig. [2 and[3}, nearly doubling the 
size of the optical disk. It is more prominent NW and SE of 
the optical disk and slightly closed to the NE, while undetected 
to the SW. The receding side of NGC 864 is narrower than the 



approaching one. The velocity field shows a symmetric pattern 
of differential rotation in the inner parts (r < 300"), while the 
kinematics of the outer parts is strikingly different. The outer 
isovelocity contours in the northern part show some twisting in 
the E direction, the last contour being nearly closed as a sign 
of a flattening rotation curve, and the southern isovelocity con- 
tours bend to the south with decreasing values characteristic of 
a declining rotation curve. The optical disk is reasonably sym- 
metric in both spiral structure and angular extent, except for a 
slight enhancement of the southern arm, and a minor extension 
of the optical disk to the SE. 

3. Modeling of the galaxy 

We have modeled the velocity-field of NGC 864 with a least- 
square algorithm based on a tilted ring model (Begeman 1987; 
the ROTCUR task in GIPSY). By this procedure we divided 
the galaxy into concentric rings, each of them with a width of 
20" along the major axis. The asymmetry of the HI distribution 
and velocity field in the outer parts lead us to model separately 
the approaching and the receding part of the galaxy. Points 
within a sector of ± 30° from the minor axis were excluded 
from the fits. The center position was fixed to the position of 
the optical center a(2000.0) = 02 /i 15 m 27: , 6 and S( 2000.0) = 
06°00'09"1 (Leon & Verdes-Montenegro 2003). The systemic 
velocity of 1561.6 km s _1 has been obtained as the central ve- 
locity of the HI spectrum at the 20% level. Expansion velocities 
were set to zero. In a first iteration the position angle, the incli- 
nation and rotation velocity have been left free, giving already 
a very satisfactory modeling of the data cube. The derived ro- 
tation curve, however, suffers from beam smoothing: since the 
beam is elongated along the major axis the integration effect 
biases the radial velocities towards values lower than the true 
values corresponding to the rotation curve. A second iteration 
has been performed, this time fixing the radial velocities at 
slightly higher values in order to correct for beam smearing, 
and leaving again the position angle and the inclination free. 
The modeled cubes for the redshifted and blueshifted parts of 
NGC 864 were combined in a single cube from which a ve- 
locity field has been obtained ( Fig. [5}- It shows a remarkable 
agreement with the observed one (Fig. [5] right), although not 
reproducing the asymmetries along the minor axis direction. 
Deriving the position angle of NGC 864 from our modeling for 
r < 300" we obtain a mean value of 23° ± 3, in very good agree- 
ment with the optical disk orientation. The inclination obtained 
as the mean between the receding and approaching sides is 43° 
± 2, which is consistent to the value given by Tully (1988) of 
45°. The modelled inclinations for both sides are almost con- 
stant and their values are always close to the mean. The major 
axis twists both in the northern and southern part of NGC 864, 
although is more pronounced to the north. There it starts at r ~ 
240" at nearly the edge of the optical disk, reaching a position 
angle of 39°. In the southern part it starts at r ~ 360" going up 
to 31°. In Fig. [5] we compare the position- velocity cut at 23° 
for the modeled (top) and observed (bottom) cube. Again the 
observations are very well reproduced. The plot also shows the 
projected rotation curve used. This figure contains the essen- 
tial peculiarities of NGC 864: asymmetric HI distribution and 
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kinematics, with a noticeable drop in rotation velocities. We 
note that the blueshifted part of the observed position-velocity 
diagram shows a faster than keplerian drop in the velocity of 
NGC 864, with a trend to reach r~° 5 values at larger radii, that 
should of course be explained by non-circular motions and/or 
projection effects. Such a phenomenon is extremely rare in HI 
envelopes. 

4. Discussion and conclusions 

Our results for NGC 864 have wider implications : 

a) The integral HI profile is symmetric in velocity, but 
asymmetric in intensity. Yet the 2D-kinematics of NGC 864 
show large scale asymmetries: there is a kinematic warp, and 
in the southern part there is an abrupt decline of the rotation 
velocity, as is evident in Fig. [6] The atomic gas here looks 
like a kinematically detached clump , evident as a secondary 
peak in the position velocity cut at radii ~ 400" to the SW. 
This region has a physical extent of around 3'x4' in the major 
axis and minor axis directions, respectively. Clearly a simple 
integrated profile analysis will conclude that the central parts 
of NGC 864 are roughly symmetric, with a value of AV20 
correlated with the luminosity of the galaxy, since it falls 
on the Tully-Fisher relation. It is the outer parts which are 
surprisingly asymmetric, but since the radial velocities in the 
approaching side are closer to the systemic velocity than those 
in the central parts, the asymmetry only manifests itself in a 
higher amplitude of the blueshifted horn of the HI profile. 

b) The outer HI envelope is large, massive, very asymmet- 
ric and presents much structure. Its HI mass is about 4.5 x 
10 9 M Q , and the HI mass associated with the steep drop of the 
rotation curve is about 2.3 x 10 8 M . No optical counterpart 
is found associated to the perturbed atomic gas. Warped HI 
envelopes are common around spirals, but few have so much 
structure and asymmetries as the one reported here. 

c) The galaxy is isolated with respect to similarly sized 
galaxies by a rather strict criterion. Yet there are 5 small com- 
panions within a projected distance of 500 kpc. The closest 
one, detected in our maps, has a dimensionless gravitational 
interaction strength (Dahari 1984) of 4 x 10~ 4 , while for 
the others this parameter ranges from 1 x 10~ 5 to 8 x 10" 7 . 
They are hence too small in mass to have caused the strong 
perturbations in the outer envelope of NGC 864. A similar 
conclusion follows from considerations about spiral forcing by 
companions (cf. Athanassoula 1984). 

We have considered several possibilities to explain the ori- 
gin of these HI asymmetries. A self-induced perturbation 
seems highly unlikely: while the pseudoring is a very standard 
feature, the outer HI ring-like structure is too large to be a bro- 
ken resonance ring, since the outer Lindblad resonance is usu- 
ally located a bit over 2 times the end of the bar (Athanassoula 
et al 1982). An external perturbation by a companion, which 
could have interacted or even been accreted, seems likely to 
have caused the SW clump and the outer ring-like structure. 



However, we can exclude an encounter with a large pericen- 
ter distance since this would necessitate a massive companion, 
which is simply not there. A small companion would have to 
come very near NGC 864 or, even better, go through it, but 
a central, or near-central, passage will (cf. Athanassoula et al. 
1997, Berentzen et al. 2003) yield results that do not match at 
all the morphology of NGC 864. The alternative left is that the 
companion crossed the equatorial plane of the target at an in- 
termediate distance, e.g. just outside the optical disc and still 
within the extended HI disc. Such a passage could have in- 
duced the warp, and if the intruder was a loosely bound gas 
rich dwarf, its gas could have contributed to the SW clump, 
which has a similar gas mass, while its stars and dark matter 
could have dispersed. To study this scenario, which is of course 
only a cartoon - and to explain the form of the ring-like den- 
sity enhancement, the mass of the excess gas and the velocity 
perturbations - full blown selfconsistent simulations, including 
both gas and stars, are necessary, but beyond the scope of the 
present paper. 
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Fig. 1. Comparison of the integrated VLA profile (solid line) and the HI profile obtained by Haynes et al. (1998) at Green Bank 
43m (dots). At the bottom we show the HI profile of the dwarf companion found in the HI map. 



Fig. 2. HI column density distribution (left) and velocity field (right) of NGC 864 and its companion, superimposed on the optical 
POSS2 red band image. The contour levels are: 5, 12, 24, 37, 49, 61, 73, 86, 98, 1 10 and 122 x 10 20 cm" 2 The velocity contours 
go from 1465 to 1660 km s in intervals of 15 km s , and are labeled each 30 km s . The beam size of 49"8 x 46"2 is shown 
in the upper left of all panels. 



Fig. 3. Grey scale map of the HI column density distribution of NGC 864. Darker regions correspond to higher column densities. 
The beam size of 49"8 x 46"2 is shown in the upper left part. The dashed line defines the approximate location of the pseudoring 
described in §|2] The insert shows the inner structure of NGC 864. 

Fig. 4. .a) HI channel plots superimposed on the optical POSS2 red band image. The beam size of 70''4 x 65''3 is shown in the 
upper left part. The contour levels are: 3, 8, 16, 32, 48, 64, 80 and 96 mJy/beam. 



Fig. 4. .b) Continued. 



Fig. 5. HI velocity modeled as explained in §0 The velocity field contours are plotted in the map in km s 1 from 1465 to 1660 
km s _1 , in intervals of 15 km s _1 . 



Fig. 6. Position velocity along p. a. = 23° for the model (top), and for the observations (bottom). The projected rotation curve is 
also plotted (dots). 



Fig. 7. Rotation curves of NGC 864 considering: a) both the receding and approaching regions, b) the approaching half (southern 
region) and c) the receding half only (northern region). The dashed line represents the rotation curve in a) without error bars. 
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